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Spectra of Galactic Cosmic Rays (GCRs) measured at the Earth are the com- 
bination of several processes: sources production and acceleration, propagation 
in the interstellar medium and propagation in the heliosphere. Inside the solar 
cavity the intensity of GCRs is reduced due to the solar modulation, the in- 
teraction which they have with the interplanetary medium. We realized a 2D 
stochastic simulation of solar modulation to reproduce CR spectra at the Earth, 
and evaluated the importance in our results of the Local Interstellar Spectrum 
(LIS) model and its agreement with data at high energy. We show a good 
agreement between our model and the data taken by AMS-01 and BESS ex- 
periments during periods with different solar activity conditions. Furthermore 
we made a prediction for the differential intensity which will be measured by 
AMS-02 experiment. 
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1. Introduction 

Models of the Heliosphere need to be very accurate in order to reproduce 
the complex structure of the solar cavity and its effect on the Cosmic Rays 
propagation. In particular simulations of CR intensity have to be compared 
with experimental values, which have been measured in the last years by 
several space experiments (e.g. AMS-01, BESS) and will be measured even 
more accurately in the near future (e.g. AMS-02 on the ISS). We have 
developed a two dimensional (radius and helio-colatitude) model of GCR 
propagation in the Heliosphere^'^ that is a function of measured values 
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of the Solar Wind velocity in the ecliptic plane (Vo) and of the neutral 
sheet tilt angle (a), as well as of stimated values of the diffusion parameter 
(Kq). This model is including curvature, gradient and current sheet drifts, 
which are depending on the charge sign of particles and magnetic field 
polarity"^. Particles modulation strictly depends on the Local Interstellar 
Spectrum (LIS), up to now not measured and supposed to be constant 
with time outside of the heliosphere. In the model we include also the 
effects of the evolving solar activity conditions experienced by a CR particle 
inside the heliosphere, due to the time spent by the solar wind to reach the 
outer border of the solar cavity. Here we deal with the differential intensity 
observed at 1 AU: we do not take into account the effects of the Earth 
magnetosphere^. 

2. Stochastic 2D Monte Carlo 

GCR transport in the Heliosphere is described by the Parker equation^'^ 

(k^M) - si7(^-^) + I'-k'MsTU) - £-S{va.)U) (1) 

where U is the CR number density per unit interval of particle kinetic en- 
ergy with time (t), T is the kinetic energy (per nuclcon), Vswi the solar wind 
velocity along the axis Xi, {v^.) is the drift velocity'', Kf^ is the symmetric 
part of the diffusion tensor and £i=(T-|-2To)/(T+To), where Tq is particle's 
rest energy. 

As demonstrated by Ito®, eq. 1 is equivalent to a set of ordinary stochas- 
tic differential equations that could be integrated with Monte Carlo tech- 
niques. The integration time step (Af) is related to the accuracy of the 
integrating process and is taken to be proportional to (r is the distance 
from the sun) to save CPU time^. We considered^° the 2D (radius and 
colatitude) approximation of eq. 1. The drift velocity is split^ in regular 
drift {vDi) and neutral sheet drift {vns)- The model is depending on some 
parameters related to the solar activity which are fine-tuned comparing 
simulations with experimental data^^. 

The perpendicular diffusion coefficient has two components, radial K±^r, 
and polar K^q. The parallel diffusion coefficient^^ is evaluated by the 
relation K\\—KQj3Kp{P){BQ/'iB). Moreover (if±)o represents the ratio 
between perpendicular and parallel diffusion coefficients: K^r={Ki_)oK\\. 
Here Ko = l-&x lO^^ cm^s-^GV-S /3 is the particle velocity, P is the CR 
particle's rigidity, the Kp{P) term takes into account the dependence on 
rigidity (in GV), ~ 5 nT is the value of heliospheric magnetic field at the 
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Earth orbit, and B is the magnitude of the Parker field-*^^. Kq is estimated 
following the procedure discussed in Section 2.2 of^^. The perpendicular 
diffusion coefficient in the polar direction K±0 has been enhanced with re- 
spect to Kj_r to reproduce the correct magnitude and rigidity dependence 
of the latitudinal cosmic ray proton and electron gradientS"'^"''. 

We use as Heliospheric Magnetic Field (HMF) the Parker field model"'^^. 
The Parker field has been modified^'' introducing a small latitudinal com- 
ponent. This modification increases the magnitude of the HMF in the polar 
regions without modifying the topology of the field. The main effect is a 
lower CR penetration along polar field lines in the inner part of the helio- 
sphere, caused by a lower magnetic drift velocity in this region, as expected 
from measured data. According to Ulysses data, for periods of low solar 
activity, we use a latitudinal dependence^^ of the solar wind velocity on the 
solar colatitude. Solar wind speed values range from Vmin — 400 km/s in 
the ecliptic plane, to Vmax — 750 km/s in the polar regions. Drift effects 
are included through analytical effective drift velocities: in the Parker spiral 
field we evaluated the three components of drift (gradient, curvature and 
neutral sheet^^ that modify the integration path inside the heliosphere. We 
adopted the approach of Potgicter & MoraaF, because, in comparison with 
other wavy neutral sheet models, it is able to reproduce the effects in both 
quite and active solar periods. 

3. Parameters and experimental constraint 
3.1. Proton LIS 

One of the inputs of our code is the proton LIS. In order to estimate the 
reliability of LIS available in literature we have considered the proton spec- 
tra measured in several periods by several experiments for rigidity greater 
than 20 GeV. In this region we can assume negligible the effects of solar 
modulation. Therefore the spectrum, which can be represented by a power 
law [N{R)r>2ogv = ^oiR/RoV] with a spectral index 7 ~ —2.7, can be 
also considered the LIS. 

We focused our analysis on four sets of data: CAPRICE- 1994^0 BESS- 
1998^^ AMS-199822and BESS-20022-^. We estimated the spectral index (7) 
and normalization constant {'^0) for each experiment, then we evaluated an 
error- weighted average on these results. We obtained: 7 = —2.77 ±0.01 and 
$0 = 15.8 ±0.6 (m^ sr MV)~^. In figure 1 we compare our results with the 
model by Burger & Potgieter""^^, where it seems to slightly overestimate the 
experimental data. We systematically repeated the above analysis fitting 
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all together the data sets obtaining: 7 = —2.77 ±0.04 and $0 = 17± 3 (m^ 
sr MV)""'^. Therefore, in the following calculations we adopt the model by 
Burger & Potgieter, corrected by a scale factor in order to fit experimental 
data. 
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Fig. 1. Fit (red line) of high energy proton spectrum measured by several experiments 
(black dots) and comparison with Burger & Potgieter^^ LIS (blue line). The spectral 
index is consistent with data, but the normalization constant is too large. 



3.2. Latitudinal Gradient 

Observations made by the Ulysses spacecraft in the inner heliosphere have 
shown that the latitudinal dependence of CR protons is significantly less 
than predicted by classical drift models^^. Our model, due to the modifica- 
tion of the i^_L6)^^ parameter in the polar regions can reproduce the gradient 
observed by Ulysses (see Fig. 3 in Heber et al.^^) between the poles and the 
ecliptic plane. We find, for the period of AMS-01 data taking (June 1998), 
a difference of 16 % between the ecliptic plane and a colatitude of ^ 30° 
from the poles. In Figure 2 a simulation of the proton differential intensity 
on the ecliptic plane is compared with the one at 30 deg from the poles 
(left panel). The latitudinal gradient of the integral intensity is shown in 
the right panel of Figure 2. 

4. Results 

4.1. Proton differential intensity 

We performed the simulations using dynamic values of Kq, a and Vsw- we 
consider present time parameters in the inner shell of the heliosphere, but 
we use the values assumed by the parameters up to 14 months back in the 
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Fig. 2. Left panel: simulation of the proton differential intensity on the ecliptic plane is 
compared with the same quantity at 30 deg from the poles. Right panel: the latitudinal 
gradient, i.e. the integral intensity for E > 200 MeV vs the solar colatitude. Simulations 
have been done using solar parameters occurred in June 1998. 



past at the heliopause. Results are shown in Figure 3. Simulated differential 
intensity with dynamic values shows a very good agreement with measured 
data, within the quoted error bars. This happens in periods with low solar 
activity and A>0, in comparison with AMS-Ol^^, BESS-QT^i and BESS-99^i; 
in periods with high solar activity, in comparison with BESS-00^^; and in 
periods with a lower solar activity and A<0, in comparison with BESS- 
02^^ and BESS-04^''. This means that our description of the Heliosphere 
improves the understanding of the complex processes occurring inside the 
solar cavity. 



4.2. AMS-02 Predictions 

Our simulation code has been used to predict CR differential intensity for 
future measurements. The assumption is that diffusion parameter, tilt an- 
gle and solar wind speed show a near-regular and almost periodic trend. 
The periodicity occurs after two consecutive 11-years solar cycles. Using 
SIDAC (Solar Influences Data Analysis Center) data we selected periods 
with a nearly similar solar activity conditions and same solar field polarity 
of the simulation time: therefore approximately 22 years in advance. We 
concentrate our simulations on the AMS-02 mission, that will be installed 
on the ISS in February 2011, for a period approaching the solar maximum: 
January 2012. Results are shown in Fig. 4. 

The Sun is currently in an unpredicted long duration solar minimum, 
that forced scientists to review all their estimations for the next solar cycle. 
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Fig. 3. Comparison of simulated galactic proton 
experimental data: AMS-Ol^^ (1998) and BESS21'24 
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therefore this prediction could be object of revision in the next future. 
AMS-02 is expected to cohect, in a fe'w years of operation, more than 10"'^° 
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Fig. 4. Prediction of modulated proton differential intensity as will be measured by 
AMS-02 experiment in January 2012. 

protons with energy > 1 GeV, and ^ 10® with energy > 1 TeV. 
5. Conclusions 

We buih a 2D stochastic Monte Carlo code for particles propagation inside 
the heliosphere. Our model takes into account drift effects and shows a good 
agreement with measured values, in periods with positive as well as with 
negative polarity. Proton spectra, as predicted by the model, are decreasing 
with increasing tilt angle and solar wind velocity. We use as input LIS 
the model published by Burger & Potgieter^^, corrected in order to fit the 
high energy measured intensity. Recent measurements have pointed out the 
needs to reach a high level of accuracy in the modulation of the differential 
intensity, in relation to the charge sign of the particles and the solar field 
polarity^®. This aspect will be even more crucial in the next generation of 
experiments. 
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